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Abstract 

In this paper we illustrate the dynamics of the instanton represen- 
tation in the description of vacuum GR in minisuperspace for unden- 
sitized variables. We uncover a new class of general solutions in both 
the degenerate and the nondegenerate sectors of the theory. Addition- 
ally, the individual sectors are preserved under Hamiltonian evolution. 
Finally, we present an algorithm for constructing general solutions by 
expansion about the isotropic sector of the instanton representation. 
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1 Introduction 



We will acquire some intuition regarding the classical dynamics of the in- 
stanton representation of Plebanski gravity, starting in this paper with a 
cursory analysis of minisuperspace dynamics in undensitized momentum 
space variables. The action for the instanton representation is equivalent 
to that in the Ashtekar variables when restricted to nondegenerate metrics. 
In this paper we will examine the dynamics of both the degenerate and the 
nondegenerate sectors of the theory. The instanton representation on the 
full phase space fi/ n st is given by the first order phase space action 

Ilnst = J dtj^ d 3 x (VaeBiA* + ^ ae B\ D t A a - N^H^j . (1) 

= (N,N l ) are respectively the lapse function and the shift vector of 
general relativity, and = (H, Hi) are the Hamiltonian and the diffeomor- 
phism constraints, given by 

H = (detBJ^Vdit^A + tr^^ 1 ); H t = e tjk BiB^ ae . (2) 

The object ^> ae is a SO(3, C) x SO(3, C) valued matrix known as the CDJ 
matrix [1], and is of mass dimension [* a e] = _ 2. B\ is the magnetic field 
for a self dual SO(3, C) gauge connection Af, 1 where 

B\ = e^djAl + l - e ^f abc A)Al. (3) 
Under the CDJ Ansatz 

vp- 1 = (a-^Bl (4) 

then (1) reduces for (deti?) / and det^ ^ to the action for GR in 
the Ashtekar variables where a l a is the densitized triad (see e.g. [2], [3], [4]). 
Hence (1) is a formulation of general relativity in which the CDJ matrix ^ ae 
is regarded as a fundamental dynamical variable. 

To study the minisuperspace dynamics of the instanton representation, 
we must now reduce (1) to minisuperspace. Minisuperspace is defined as 

x Our notation is that symbols from the beginning of the Latin alphabet a, b, c, ... 
denote internal 5*0(3, C) indices, while from the middle of the alphabet i, j, k, . . . denote 
spatial indices. 
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the sector of the full theory where all variables are spatially homogeneous. 
This means that all spatial gradients must be set to zero, which is unlike 
the usual definition of minisuperspace which uses Bianchi symmetry groups 
[5]. Hence in minisuperspace as we have defined it, (3) reduces to 

B\ = (detA^A- 1 )^ detB = (detA) 2 . (5) 

We have used in (5) the fact that the structure constants f a b c for SO(3, C) 
are numerically the same as the three dimensional epsilon symbol e a b c , in 
writing the determinant. We must now reduce the constraints (2) to min- 
isuperspace. The Hamiltonian constraint H is given by 

H = (detA)Vdettf(A + tr^ 1 ) (6) 
and the diffeomorphism constraint Hi is given by 

Hi = (detA) 2 (detA)' 1 Afy d = (detA)^d, (7) 

where ipd = edbf^bf is derived from the antisymmetric part of ^fbf- A direct 
way to obtain the Gauss' Law constraint for minisuperspace is to obtain 

V ae BlD t A a = VaeBl&AZ + f abc A^A c ) = * ae (det A) (A^fj^A^ (8) 
and then vary (8) with respect to Ag, yielding 

G a = (detA)Vv (9) 

In minisuperspace the both the Gauss' law and the diffeomorphism con- 
straints depend linearly on ip^, and are therefore redundant. 

The last remaining object needed is the canonical structure which deter- 
mines the canonical one form, which is given by 

X ae = B l e Af = (detA)(A~ 1 ) i e Af . (10) 

While the velocity X ae is defined by (10), the issue of the existence of X ae 
as a global coordinate on configuration space Tj nst arises. 2 As we will show, 
the existence or nonexistence of X ae is not relevant as far as the dynamics 
are concerned, since we will be able to formulate the Hamilton's equations 
of motion using only the velocities X ae without actually making use of X ae . 

2 The Soo one forms SX ae were first introduced in [6] and [7]. 
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The organization of this paper is as follows. Section 2 treats the nonde- 
generate case, computing the general solution to the equations of motion, 
and section 3 derives the induced time development of the spacetime met- 
ric. Sections 4 and 5 recompute the equations of motion for the degenerate 
case, and construct some new solutions. Section 6 introduces a method for 
constructing a general solution by expansion about the isotropic sector. 
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2 Setting the stage with nondegenerate vacuum 
general relativity (X = detA ^ 0) 

The action for general relativity in the instanton representation in minisu- 
perspace with Lorentzian signature is given by 



nst 



+(detA) [iivVdet* (A + tr^- 1 ) + (WAf - ^)/ d6/ * 6/ ] ) , (11) 
The Hamilton's equations of motion for X ae are given by 



X° 



BlAf 



5H 



-Af(detA)v / det*(^- 1 ^- 1 ) ae + (deU)(JV^ - 9 d )f dae ; 
and for the CDJ matrix \Ev by 



(12) 



A TT r- 

= -gx^ = ~ [iSaeidetA^HiN] 
+8 ae (N*Af - 6 d )f dbf * bf 



(13) 



Since the Cauchy development of GR should be consistent with the initial 
value constraints, then the constraints can be applied wherever they appear 
in the equations of motion. For det£> = (det^4) 2 / the only constraint 
contained in (12) is the Hamiltonian constraint H ~ 0. 3 On the other hand, 
all terms of (13) are directly proportional to constraints. 

The diffeomorphism constraint Hi = implies that fdbf^bf = 0, or 
= *(&/) is symmetric in bf. Application of this and the Hamiltonian 
constraint H = in (13) yields 



*ae = -> *oe(*) = \ae) = COUSt. Va, C 



(14) 



However, for (12), only the Hamiltonian constraint can be used. The solution 
to the Hamiltonian constraint is given by 



A1A2 



AAiA 2 + Ai + A 2 



(15) 



3 We have used the identity B l a = (detA)(A 1 )*, as well as X = trX ae = detA, for 
anisotropic minisuperspace. 
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To obtain some physical insight into the configuration space dynamics, 
let us take the symmetric part of (12), which has the same effect as making 
a gauge choice iV* = 6 d = 0. Then (12) and (13) reduce to 



x (ae) = _ iN (^ et A)r] ae - V ae = 0, (16) 

where we have defined an internal SU(2)- (8> <S77(2)_ metric rf e (not to be 
confused with the Minkoski metric rfi which has spatial indices by 



7] ae = Vdet^i^- 1 ^- 1 )^. (17) 

Equation (16) resembles the equations of motion for a free particle in classi- 
cal mechanics, travelling ostensibly through a nine dimensional configuration 
space. However, the space is actually five dimensional when restricted to 
the constraint surface defined by the initial value constraints. 

Taking the trace of (16) and dividing through by X = trX ae = detA 
since I/O due to nondegeneracy, we obtain 4 



X d\nX , . 

x- — = ~"> N (18) 

where we have defined rj = 5 ae r] ae as the trace. Since X is globally a holo- 
nomic coordinate on Tj nst , then equation (18) directly integrates to 



X(t) = (detA(t)) = X e- ir > N ^ dt ' (19) 

where Xq = X(0). Note the mass dimension [77] = 1 which cancels the 
negative mass dimension of time [t] = —1 so that the argument of the 
exponential is dimensionless. Substituting (19) back into (16) we obtain the 
equation for the velocity X ae 



X(«0 = -iNr, ae (x e- ir 'fi N W ( #) = — ^ (x e _ir? N (t')dt'^ 

(20) 

The right hand side of (20) is a total time derivative since according to (14) 
and (17), r] ae is constant in time. This implies that for the equation to 
make sense the left hand side must as well be a total derivative. However, 
there is no coordinate X ae on configuration space, since 5X ae = B l e 5Af <j£_ 
/\ a (Ti nst ) is not an exact one form except for its trace. However, there exist 

4 Note that X is actually the Chern-Simons functional for minisuperspace, identified 
in [7] as the candidate for a time variable on configuration space in quantum cosmology. 
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configurations for which a densitized X ae can in some sense be defined. For 
instance, for diagonal connections Af = 5? A® one has 

X = (detA) = A\A\Al- B\ = A 2 2 A 3 3 ; B\ = A\A\; B| = A\A 2 2 (21) 
and (16) reduces to 

A\A\A\ = -iNr, n (A\Jl*Al) — ► = -i*V U , (22) 
and likewise for A\ and A\. Then (22) integrates to 

for / = 1,2,3. The result is that all components of the connection evolve 
with respect to X, seen as a time variable on configuration space. 

2.1 Dynamics of the spacetime metric 

Since the instanton representation is a metric-free description of gravity, then 
the spacetime metric is a derived quantity. The lapse function N(t), 
can still be chosen arbitrarily, but the spatial 3-metric h^ is determined 
dynamically through the evolution of the instanton representation phase 
space. The spatial three metric is given by 



det# VdetB/ deter' V ' 

The relation to rj ae in (17) stems more directly from the covariant form 



hij = (det*)(^- 1 ^- 1 )^(B- 1 X e (B- 1 )J(detB), (25) 
whereupon the following relation can be written 



rf e = Vd^^- 1 ^- 1 )^ = (dettt)" 1 / 2 ^ ■ f^f ) . (26) 

V deti? / 

Making use of the minisuperspace relations B % a = (&etA){A~ l ) % a and det-B = 
(detA) 2 = X 2 , and v / det* = (det??) -1 , then (25) can be written as 



hij ~ hij(t; A) = (det?7)-y e (A)^(t)^). (27) 
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Observe that h^ has acquired the label of rf e G GL(5,C), and also 
depends on Af. 5 The part of hij which can always be unambigously specified 
is h = det(hij) = (detrj)~ 2 X 2 , whose time evolution is given by 

y/htfj = (detr ? )- 1 X(0)e- i ^o JV (*')<it'. ( 28 ) 

Hence in the general solution for the nondegenerate case even h, acquires 
the label of X ae £ GL(5, C). One could then construct a spacetime metric 

ds 2 = -N 2 dt 2 + (detr?)- V e cj a <g> u e , (29) 

whence the Ashtekar potential becomes absorbed into the definition of 
the one forms uj 1 , given by 

u a = A^(dx i + N i dt). (30) 

While (30) resembles the invariant one forms defined for Bianchi groups [5], 
this is not the case since A® does not satisfy the Maurer-Cartan equation 
since it is not a flat connection. 

For the configuration (23) one obtains the time evolution for the spatial 
3-metric 

/ X (t\ \ 2-nf-t In 

hijit) = SifSjfidetvrWHA^ij^) • ( 31 ) 

which also evolves in relation to X(t), seen as a time variable on config- 
uration space. This can be written explicitly in terms of metric variables 

as 

h ij (t) = 5 if S jf h ff (0)(^y f,/ \ (32) 

where now the determinant of the 3-metric plays the role of the time variable, 
and the components of the metric evolve with respect to it. From (32) one 
finds the physical interpretation of the components of the connection Aj in 
terms of hij , which fixes its value at t = to 

h ff (0) = (det^-y^O)) 2 . (33) 

5 According to the previous section, the information encoding the evolution of specific 
components of Af, is tied up in the combination SX ae = BISA". The trace of X ae , 
namely the Chern-Simons invariant, undergoes a well-defined evolution (19) from which 
some information regarding A^(t) can be inferred. 
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To obtain a real section of GR, one imposes reality conditions on the in- 
stanton representation variables so that the 3-metric is real. For example, 
we must have ty* /rj real for each /, namely that 

Im(r] u r] 22 ) = Im^rf 3 ) = Im(r] 33 r] U ) = 0. (34) 

The real parts of 77 11 , rj 22 and rj 33 and one imaginary part, say Im(r] 33 ) are 
freely specifiable, which fixes /m(ry 11 ) and Im(r] 22 ). In order for the metric 
to be real at t = 0, then (33) in turn requires that each component of the 
connection be either pure real or pure imaginary. 
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3 The degenerate sector (detA=0) 



Having obtained a solution for the nondegenerate sector, let us now examine 
the degenerate sector (detA) = 0. This implies that the metric is also 
degenerate, since h = (det^4) 2 (det^'). While this may be the case, as we will 
show, the time evolution of the CDJ matrix is still well-defined. Recall the 
expression for the Hamiltonian constraint 



H = (detA)Vdet*(A + trl^ 1 ) , (35) 

which includes detA as part of its definition, therefore when detA = (35) 
is trivially satisfied with no restrictions on ^ ae . This means that the CDJ 
matrix fy ae is now free to evolve in time. Let us now revisit the equations of 
motion (12) and (13) under the condition of degeneracy. Starting with the 
equation for X ae , requoting (12) for completeness, 



X ae = BtA? 
-iV(detA)Vdet1'(^- 1 ^- 1 ) c 



(deU)(JVM?-0/dae- 



(36) 



Since detA = 0, then (36) the trace of (36) vanishes. This implies that 
the trace X is numerically constant, and since X(t) = detA = = Xq, it 
remains degenerate for all times. Therefore, the degenerate case remains a 
distinct sector of vacuum GR which cannot be bridged. 6 

Moving on to the equation of motion for * ae , we have from (13) that 



<H ae = -i6 ae NVde~M{A + tr*- 1 ) - MiVM? " ^fdbf^bf- (37) 

The Gauss' law and the diffeomorphism constraints imply that = 

is symmetric in b, /, therefore the second term of (37) vanishes and we are 

left with the equations 



4> ae = -i5 ae NVdet^(A + tr^ 1 ) ; X ae = 0. (38) 

Comparison of (38), corresponding to (detA) = with the analogous equa- 
tions (16) for (detA) / reveals the following contrast. Whereas for 
(detA) 7^ the momentum variable ^> ae was constant in time while X ae 
was nontrivial, we see for (detA) = that it is X ae which vanishes while 

6 Hence, the configuration space variables cannot evolve in time from the X — sector 
into the X ^ sector, since the equations of motion are assumed to hold for all time. 
Consequently, topology change is precluded, unlike in [11] and works by other authors. 
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\^ ae inherits a nontrivial time evolution. The 'configuration' and the 'mo- 
mentum' space have essentially 'exchanged' roles as a consequence of the 
degeneracy condition. 

Equation (38) states that the time derivative of \l/ ae is an isotropic ma- 
trix, which means that the off-diagonal parts are numerical constants. Since 
^?[ ae ] = on account of the kinematic constraints, we may assume that ^ ae 
is a symmetric matrix of the form 

/ a(t) W V \ 

*ae(t) = W b(t) U 

V V U c(t) J 

where U, V and W are arbitrary numerical constants. Without loss of 
generality we can set U = V = W = 0, since ^> a e can always be diagonalized 
if it is nondegenerate. Note that these initial conditions are preserved on 
account of the equations of motion, and the nontrivial equations of motion 
then reduce to 

a = b = c = iNVabc(A + - + - + -). (39) 
V a b c) 

The time derivates of the diagonal elements are equal, which implies that 
these elements must be equal, within numerical constants, to each other 

b(t) = a(t) + k 1 ; c(t) = a(t) + k 2 (40) 

for arbitrary constants k± and &2 . In the general case the equation of motion 
can be integrated 

ra(t) / / 1 1 1 \ 

/ ( a(a + kl )(a + k 2 )) ~ 1/2 ( A + - + — tt + —rr) 
Ja V a a + ki a + k 2 J 

The left hand side of (41) can be written in closed form in terms of known 
functions, but we do not display the result it here. Let us rather use a 
compact notation to decribe the integral 

I A {a;k 1 ,k 2 ) = i f N{t')dt' = r (t). (42) 
J o 

One can then in principle invert (42) to find the evolution of a = a(r) as a 
function of the 'time' r, which depends on the choice of lapse function N 
and is labelled by the constants k\ and k 2 . 



f N(t')dt'. (41) 
Jo 
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3.1 A few simple cases within the diagonal sector 

Let us now illustrate a few simple examples for the degenerate case for which 
short expressions can be written. 

Case (i): Isotropic case with A / 0. In this case we have k\ = £;2 = 0. The 
equation of motion reduces to 



a = iNaV 2 (k + -^ 
Equation (43) directly integrates to 



(43) 



lA(a;0, 0) = r(i); a(r) = —tan 



tan 



l ( /Aa y/3A 



( 



For A = we have 



3 \ 2 



o(r) = ( v / «o+ 2 T ) 



3 



(44) 



(45) 



Case (ii): One diagonal degree of freedom. In this case k-i = 0, k\ = k where 
k is an arbitrary numerical constant. The we have 



(2 1 \ 
ah — i ) 
a a + kJ 

Equation (46) integrates to 



(46) 



l\(a;k,0) 



..tanh" 1 (3^3) -r+tanh-^^S) 
kA - 3 - r 2 



r(47) 



where we have defined the dimenionless constants 



= \JkA-3- Vfc 2 A 2 - 2fcA +^5 r+ = V^A - 3 + ^k 2 A 2 - 2fcA+^. 



(48) 



While solved in closed form, (49) is nontrivial to invert for a = a(t), but 
nevertheless the solution is implicit. The A = case leads to the relation 



I (a; k, 0) = - Va + k - J -tanh" 1 Us (l + ^) ) = r. 



(49) 



Likewise, the relation (49) is implicit but still nevertheless illustrates the 
integrability of the system. 
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Case (iii): General diagonal case for A = 0. This is given by 



T( l, y\ f a(t) V<a + ki)(a + k 2 ) 

I {a;k 1 ,k 2 ) = 2 . . da = r 50 

7 ao 3a 2 + 2a(fci + fc 2 ) + kik 2 

The relation (50) can as well be integrated in closed form in terms of known 
functions, though we do not display the final expression here. 
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4 The degenerate general case 



We now treat the general solution for the degenerate case, solving the equa- 
tion of motion 



4> ae = -id ae NVdeW (A + tr^ 1 ) ; X ae = (51) 
for a more general form of ^ ae . This is given by 

/ a W V 

*ae=\ W b U 

\ V U c 

It will be convenient to parametrize ^> ae by its isotropic and its non-isotropic 
parts, as in 



^ae = &ae<P + ^ae- (52) 

In (52), ip is the isotropic part and a ae is defined as the 'CD J deviation 
matrix'. Substitution of (52) into (51) yields 

e ae = 0; e ae = const. Va, e (53) 

The deviation matrix e ae is in general an 8 by 8 matrix of arbitrary complex 
numerical constants. 

Let us now re-examine the equation of motion for ^ ae . This can be 
written as one equation labelled by the constants e ae G GL(8,C). 

<p = -iN^detidaeip + e ae ) (A + ti(5 ae <p + e^)" 1 ) . (54) 
To solve the equation of motion (54) we will need the determinant, given by 

det^ ae = det(5 ae (p + e ae ) = (p 3 + (p 2 e + -ipVare + dete (55) 
where we have defined 



e = 5 ae e ae = tr(e oe ); Vare = {5 bf 5 cg - S cf S bg )e bf e cg = (tre) 2 - tre 2 . (56) 

We will also need the inverse of (52). Let us assume for simplicity that the 
isotropic part ip is the largest component. 7 Hence, ip > e ae Va, e, with ip ^ 0. 
Then we have 

7 The trace is the dynamical component and will vary in time. Therefore, the expansion 
is good only for t = r such that tp(r) > max{e ae }. Since the non-isotropic components 
are numerical constants, then V t such that ip(t) > e a e Va, e do not hold, then a singularity 
exists and one can must attempt to find the solution for different e ae . Hence we have 
shown that singularities can develop only when considering degenerate configurations. 
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®ae = (<W + <^ae)~ l = f^Saao (~ l T ^ e a l)ai £aia 2 ■ ■ ■ ^a n - ia „) $a n e (57) 

The trace of (57) is given by 



n=0 



trtt 1 = ( p x ^(-l)> n e aai e aiQ2 ...e a „_ ia 



(58) 



n=0 



Likewise, the square root of the determinant can be expanded as in 



n — t- inf. e V are deteN 1 /- 
Then one can use the infinite binomial series expansion 



(59) 



e Vare detexV 2 ^ (-l) n (2ra)! re Vare dete\™ , . 
+ — + ^ =g (l-L)nl^ (^ + + — ^ - (60) 



We must now put (58) and (60) back into (54). 



4.1 Asymptotic expansion about the isotropic sector 

We are now ready to write down the general solution for the degenerate case 
to all orders. The equation of motion for ip can be written as 



(A + -)v?~ 3/2 vi = -iN'Efo, e) 



where we have defined a 'correction' factor for non-isotropy E 



E(<p,e) 



h(<P,eae), 



where 



(61) 



(62) 



n=l 



( e Vare detex 1 / 2 



(63) 
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The left hand side of (61) reduces to the isotropic case considered in (43). 
Upon integration, we have 

J*(A + ^<p-V 2 dy> = -i J* N(t')E(<p(t'),e)dt' (64) 
which can be written in the form 

Equation (65) further simplifies to 



We can now generate an asymptotic expansion about the exact isotropic 
solution by fixed point iteration procedure, assuming that the right hand 
side of (64) is small. First, choose a particular state e ae , thought of as 
a constant vector in an eight complex dimensional space. Next, define a 
sequence of functions of time (f n (t), such that ip = tpo, which is the initial 
value of the isotropic part of the CD J matrix at time t = 0. 8 Define sequence 
(p n , such that 



¥>n+l(*) 



yftanftan- 1 (^) + ^ jT Efatf), e)dt'] . (67) 



The full solution, assuming convergence of the sequence (67), is given by 
<p{t) = lim Tt _s. 00 tp n (t). Hence, if it is indeed the case that there was a phase 
of the universe in which degenerate metrics exist, (i) The universe cannot 
transition into a nondegenerate state, (ii) The existence of such a configu- 
ration is labelled by eight arbitrary complex constants for the e ae , each of 
which defines a state of the universe, (iii) The convergence of the sequence 
(67) is a necessary, but not sufficient condition for the well-definedness of 
this configuration. 9 This in turn depends upon the initial value of ipo in 
relation to the constants e ae . The theorist has free choice of e ae , within the 
range of parameters yielding a fixed point. However, <po is determined by 
the initial conditions of the universe, which can only be extrapolated from 
the cosmological data in existence today. 



This time labels the beginning of the universe, and ipo is an essential parameter of the 
theory. 

9 There is clearly a large category of cosmological evidence that dictates that the space- 
time metric today is nondegenerate. 
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5 Conclusion 



We have illustrated some of the dynamics of gravity in the instanton rep- 
resentation in minisuperspace for undensitized CD J matrix ^ ae , seen as a 
momentum space variable. Starting from this representation, which is de- 
fined independently of the existence of a spacetime metric, we have derived 
solutions which include a class of inflating spacetimes labelled by the eigen- 
values of ^>ae- We have also demonstrated that the big bang singularity can 
be avoided in this representation independently of any quantum corrections 
by applying the equations of motion to arbitrarily small times, including 
the beginning of the universe when quantum effects are presumably domi- 
nant. We have also illustrated the dynamics of the theory in the degenerate 
sector, displaying an array of new solutions. Another result of the intial 
value constraints development in these variables is that an initially degener- 
ate solution will remain degenerate for all times, which eliminates topology 
changing configurations classically. 

The present paper has demonstrated the anisotropic minisuperspace sec- 
tor of the model. We have also provided a prescription for obtaining general 
solutions by asymptotic expansion about an isotropic solution. It is hoped 
that this algorithm can be of use as well in numerical treatments of general 
relatity. While the configuration space variables X ae which would ordinarily 
be conjugate to ^> ae are not globally integrable for generic configurations, 
there is certain information which can be inferred from specific combina- 
tions of these variables which corresponds to a well-defined time evolution. 
A future paper along these lines will investigate the structure of X ae , using 
a version of the instanton representation where the variables are densitized, 
and as well carry out the analogous computations of the present paper. 
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